NACA TN No. 1630 / 


'Wo t 

i 


rJ^J^ 


Uj 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 

No. 1630 

A GENERALIZED THEORETICAL INVESTIGATION OF THE 
HYDRODYNAMIC PITCHING MOMENTS EXPERIENCED 
BY V-BOTTOM SEAPLANES DURING STEP-LANDING 
IMPACTS AND COMPARISONS WITH EXPERIMENT 

By Benjamin. Milwitzky library COPY 

Langley Memorial Aeronautical Laboratory c 

Langley Field, Va. JUN 1§ 

lanolwfSovhw'N'a 


MB DOCUMENT ON LOAN FROM THE FU IS Of 


NATIONAL ADVISORY COMMITTEE FOR AER0NA"T!CS 
LANGLEY MEMORIAL AEROM’|T<-. V LA Yl'fSJORT 
LANGLEY FIELD, HAMPTON, Vs*f I*(A 

fctETURN TO THE ABO'/E AOORE'l* 




Washington 
June 1948 


REQUESTS FOR PU*3L"7A f O'-* S'* 1 Bfc -tO 

AS FOLLOWS: 


JH1S DOCUMENT ON LOAN FROM THE FILES OF 

NATIONAL ADVISORY COMMITTEE FOR AERONA'TIC 
LANGLEY MEMORIAL AERONAUTICAL LABORATORY 
LANGLEY FIELD, HAMPTON, VIRGINIA 

RETURN TO THE ABOVE ADDRESS. 


NATIONAL ID 
17 ’i S* : 
WASHINCTTV 


,t f *+ .tty vi TjCS 


C- 


REQUESTS FOR POBUCAlTONS SHOULD 8C AOURCSTCO 
AS FOLLOWS: 

\ national advisory committee for aeronaut^ 

‘ 1714 F STREET, N.W., 

WASHiNGTOfl 25. D.C. 




NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO. 1630 

A GENERALIZED THEORETICAL INVESTIGATION OF THE HYDRODYNAMIC 
PITCHING- MOMENTS EXPERIENCED BY V-BOTTOM SEAPLANES 
DURING &ETT— LANDING IMPACTS AND 
COMPARISONS WITH EXPERIMENT 
By Benjamin Milwitzky 

SUMMARY 


A theoretical Investigation Is made of the hydrodynamic pitching 
moments experienced "by V-bottom seaplanes during step— landing impacts. 

The entire immersion process is analyzed from the instant of initial 
contact until the seaplane redounds from the water surface. In the 
analysis the primary flow about the immersed portion of a keeled float 
or hull is considered to occur in transverse flow planes and the virtual- 
mass concept is applied to calculate the reaction of the water to the 
motion of the seaplane. 

t 

It is shown that the pitching moment and the center-K3f— pressure 
location may "be represented in generalized form by means of dimensionless 
variables which take into account th.9 effects of such factors as the 
seaplane weight, angle of dead rise, trim angle, and initial velocity. 

The variation of these coefficients during an impact ‘is governed by the 
magnitude of a single parameter, called the approach parameter k* that 
is determined by the trim and the initial flight— path angle and permits 
reduction of the impact variables to a common basis. 

Equations are presented from which the variation of the pitching- 
moment and centeiwDf— pressure coefficients have been calculated for a 
wide range of conditions extending from impacts along shallow flight 
paths approaching planing to very steep impacts where the resultant 
velocity is normal to the keel. Solutions are also presented for the 
conditions which exist at the instants of maximum acceleration, maximum 
moment, maximum draft, and at the instant of zero draft during rebound. 
The results show that the maximum hydrodynamic pitching moment is 
attained slightly after the maximum acceleration is reached and occurs 
prior to the attainment of the maximum, draft. 

The analysis also shows that the pitching moment about the st9p at 
any given stage of the impact process is independent of the angle of 
dead rise and that the center of pressure in an impact is located at a 
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distance only slightly greater than one— -third the vetted length forward 
of the step. The most forward center— of— pressure locations, in terras of 
the vetted length, are attained at low values of- the approach parameter ic 
(■which, for a given trim angle, correspond to steep flight paths) hut 
the centeiwxf— pressure distance converges to a limiting value of exactly 
one— third the vetted length forward of the step as the planing condition 
is approached. 

Comparisons of theoretical pitching-moment time histories and maximum 
values of the pitching moment with expert mental data obtained in the 
Langley impact basin with floats of 30° and 40° angles of dead rise 
indicate that the calculated results are in good agreement with the 
measured values. 


INTRODUCTION 


Because of the evident need for a more rational foundation upon 
which to base water-loading requirements for the design of modern sea- 
planes, an extensive theoretical and experimental research program has 
been undertaken in the field of hydrodynamic impact loads. The theoret- 
ical phase of the program has resulted in references 1 and 2, which 
deal with the analysis of impact loads for V— ' bottom seaplanes. Refer- 
ence 1 presents a critical survey of previously published impact theories, 
which are shown to be applicable only to impacts where the resultant" 
velocity is normal to the keel, and in addition considers the effect of 
the velocity ccmponent~parallel to the keel. Reference 2 provides a 
generalized analysis, with experimental confirmation, of the motions and 
overall hydrodynamic loads experienced by V— bottom seaplanes throughout 
the course of step— landing impacts and shows that the motion and time 
characteristics of such impacts may be represented by means of dimension- 
less variables which take into account such factors as the seaplane 
weight, dead— rise angle, trim angle, and initial velocity. The variation 
of these coefficients during an impact was shown to be governed solely 
by the magnitude of a single parameter, called the approach parameter k, 
•which is determined by the trim and the flight— path angle at the instant 
of initial contact with the water surface. 

The present paper represents a continuation of the work of refer- 
ence 2 and applies the theoretical coefficients and equations introduced 
by that paper to the development of the nondimens ional hydrodynamic 
pitching-moment and center— of— pres stir e coefficients for step— landing 
impacts. Tbe entire immersion process is analyzed from the instant of 
initial contact until the seaplane rebounds from the water surface. In 
the analysis the primary flow about the immersed portion of a keeled 
seaplane float is considered to occur in transverse flow planes and the 
virtual-mass concept is applied to calculate the reaction of the water 
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to the motion of the seaplane. The results of the analysis are presented 
in the form of charts which may he used to determine the hydrodynamic 
pitching moments and the center— of— pressure location on a seaplane at any 
time during an impact as veil as the pitching moments, the center*-of— 
pressure location, the state of motion of the seaplane, and the time after 
contact at the particular instants of maximum pitching moment, maximum 
acceleration, maximum draft, and at the instant of zero draft during 
rehound. The applicability of the theoretical results is illustrated hy 
comparisons with experimental data obtained in the Langley impact basin 
with a Y-^bottom float (fig. ,l) having an angle of dead rise of 30° at the 
step and a similar float of 40° angle of dead rise. 


SYMBOLS 


A 

a 

b 

c 

F 

g 

l 

M 

n iv 

P 

8 

t 

Y 

W 

X 


hydrodynamic aspect ratio 

longitudinal distance between step and center of moments 

longitudinal distance between center of gravity and center 
of moments 

distance, measured normal to longitudinal axis, between center 
of gravity and center of moments 

force 

acceleration due to gravity 
wetted keel length 
pitching moment 
two-dimensional virtual mass 

impact load factor, measured normal to water surface 

longitudinal distance between step and center of pressure 

distance from given flow plane to foremost immersed station 
along keel 

time after contact 

velocity of seaplane 

weight of seaplane or total testing weight 
distance parallel to water surface 
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draft of keel at step, normal to water surface 
penetration of a given station, normal to keel 
angle of dead rise 

flight-path angle relative to water surface 
mass density of water 
trim angle 

function governing variation of virtual, mass with dead rise 
end-flow correction to total hydrodynamic load 
end- flew correction to hydrodynamic pitching moment 
Sub scripts: 
f float 

h horizontal or hydrodynamic 

rnn.-y maximum 

o initial conditions or referring to point o 

r resultant 

s at the step 

T total 

v vertical 


k 

y 

z 

0 

7 

P 

T 

f(P) 

0(A) 

0l(a) 


Dimensionless Variables 


Pitching-moment coefficient based on vertical velocity 
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Pitching-moment coefficient "based on resultant "velocity 


„ . - 1b /g\ 

W ^(A) 


Center— of— pres stir e coefficient 


0(X) 

°cp = P sin t — — 

0l(A) 


r, 


jf(p)J 2 0( A )P rt 


6 sin r cos^r 


kl/3 


Batio of center— of -pres sure distance to vetted length 


p 0(A) 
r " * t(A) 


Load— factor coefficient 


° V e 


r Q l\l/3 

6 sin t cos^r [\ 

(f(p)] 2 0(A)pjt j^/ 


° J = 15 1 


Vert leal—" velocity ratio 


To 


Draft coefficient 


° 4 - r If A F ^- — ^ }) 

V (6 sin t cos t J 


[f (p)] 2 0(A)prt| 
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Time coefficient 


°t - tfo (f 


if (p)] g 0(A)pjt 
1 6 sin t cos 2 t| 


1V3 


Approach parameter 


* = t£t ooe (T + rJ 


ANAIZSIS 

Basis of Analysis 


A detailed discussion of the physical concepts upon -which the 
analysis is "based is given in references 1 and 2. Briefly, the flow 
about a slender Immersing shape such as a keeled seaplane float or hull 
is assumed to occur in transverse flow planes which may he considered 
fixed in space and oriented normal to the keel. (See fig. 2.) Because 
of the absence of a satisfactory three-dimensional theory, the motion 
of the fluid in each flow plane is treated as a two-dimensional phenomenon. 
In order to account for the effects of end flow that exist in the three- 
dimensional case, the total force on the seaplane, which is obtained by 
integrating the reactions of the individual flow planes In contact with 
the hull, is reduced by the application of an aspect— ratio type of 
correction. 

In any flow plane the momentum imparted to the water Is determined 
solely by the growth of the float cross-sectional shape Intersected by 
the plane and may be expressed as the product of the virtual mass 
associated with the immersed cross section and the velocity of penetration 
into the plane. After the step has passed through a given flow plane the 
Intersected cross section ceases to exist and the plane becomes part of 
the downwash where it remains thereafter, independent of the subsequent 
progress of the impact. . 

In potential flow the two-dimensional virtual, mass of any float 
cross section is determined by the shape immersed in the flow plane. In 
the case of V— shaped, cross sections. If the chines are not immersed, the 
flow patterns at all degrees of penetration are models of each other 
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and the -virtual mass is proportional to the square of the penetration. 
The constant of proportionality is determined "by the dead— rise angle p. 
Thus the two-dimensional virtual mass of any V— shaped cross section may 
he expressed hy 


XV = |f(f3 )] 2 ^ z 2 


( 1 ) 


-where the quantity f ( 3 )z represents the radius of an equivalent semi- 
cylinder of water constituting the virtual mass of the immersing shape. 

On a two-dimensional basis , the momentum of the water contained 
within any flow plane is 


V = 




( 2 ) 


The reaction of the water within the plane is therefore given hy 

F = jf(3)J 2 §£■ (z 2 z + 2 zz 2 ) ( 3 ) 

The total force on the seaplane, which acts normal to the heel, is 
obtained hy integrating the forces contributed hy the individual flow 
planes in contact with the hull. Since conventional floats and hulls 
are essentially prismatic for an appreciable distance forward of the 
step, the integration is performed under the assumption that the immersed 
portion of the hull has constant cross section and that the trim remains 
constant during the short duration of the impact j this integration gives 


F t = 


[f (3)] 2 0(a) 


prt 

2 tan r 



w 


where the quantity 0(A) is applied to the calculated total force as a 
correction for the effects of end flow, -which exist in the three-dimensional 
case and is determined by the geometry of the immersed portion of the 
seaplane . 

In reference 2 an analysis of the equations of motion resulting 
from the application of equation (^) to the condition where the wing lift 
is equal to the weight of the seaplane showed that the motion and time 
characteristics of such impacts xnay be represented in generalized form 
by means of the following dimensionless variables: 
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Load— factor coefficient 


n l w .S (W }& sin T cos 2 t I 

y G 2 U ) [f(p)]^(A)p* 


l \!/3 


(5) 


where 


Draft coefficient 


ni w ~ g 


C d = y 


£ f[ f (PQ 2 0( A ) pTC 

^ 1 6 sin t cos^r 


-/3 


( 6 ) 


Time coefficient 


„ 1 e J &(P)] 0 U)p* 

* " 70 1~Bln“';'ca2r 


VL /3 


(7) 


Yertical-veloclty ratio 


( 8 ) 


'.o 


It was also shown that the variation of these nondlmensional 
quantities during the course of an impact is governed solely by the 
magnitude of the approach parameter : 


K 


sin t 
B ln ~0 


cos (r + 


^o) 


(9) 


which may be considered a criterion of impact similarity and depends only 
on the trim and the flight— path angle at the instant -of initial contact. 
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For a given value of k, therefore, the respective variations of the 
preceding dimensionless variables during an impact may each he represented 
by a single curve, regardless of -what the dead— rise angle or weight of 
the seaplane, the attitude, or initial velocity may be. Similarly, there 
is a single variation with k of each of the dimensionless variables 
representing the state of motion and the time at any given stage of the 
impact. Figure 3 is a graph of equation ( 9 ) and shows the variation of k 
with trim and flight— path angle. 

In the present analysis the theoretical coefficients and equations 
presented in reference 2 are applied in the development of the non- 
dimensional pitching-moment and center— of— pressure coefficients. As a 
result of this generalized treatment, the number of independent variables 
and the number of cases for which solutions are required are considerably 
reduced while the presentation and correlation of data for the entire 
range of seaplane and flight parameters is greatly simplified. 

Although the application of the theoretical results to impacts of 
particular seaplanes requires the definition of f(|3) and 0(A), the 
generalized solutions of the motion and time characteristics of step 
landings in terms of the foregoing dimensionless variables are valid 
regardless of how these functions may be defined. The function f(p), 
which governs the variation of virtual mass with dead rise, has been 
evaluated from the results of an iterative solution for the force on a 
two-dimensional V— shape immersing with constant velocity, reported by 
Wagner in reference 3* This variation is given by 


f (P) = 20 ~ 1 


( 10 ) 


The end— loss correction 0(A) may be approximated from the results of 
tests conducted by Pabst with vibrating plates in water (reference 4). 
Application of these results to the V— bottom seaplane gives 


0(A) 


tan t 
2 tan 3 


(H) 


Although the validity of the relationships given by equations (10) 
and (ll) has not been experl mentally verified for very low dead— rise 
angles, the application of these functions has yielded calculated 
results that closely agree with extensive test data obtained in the 

n 0 

Langley impact basin with floats of 22j!j- , 30°j and 40° angles of dead 
rise and reported in references 1, 2, and 5- 
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Pitching Moments 

As shown by equation ( 3 ), the force contributed by a given flow- 
plane may be considered to arise from two sources: namely, the inertia 
reaction to the acceleration of the virtual mass and the additional 
change in momentum accompanying the expansion of the virtual mass with 
penetration into the flow planes. Under the assumption of two-dimensional 
flow within the flow planes, for a prismatic hull at positive trim, the 
inertia reaction follows a quadratic variation along the keel while the 
force due to the expansion of the virtual mass is linearly distributed. 
(See fig. 2.) The shape of the longitudinal distribution of the total 
force during an impact is, of course, determined by the relative 
magnitudes of the component distributions. For steady— state planing, 
however, since there is no acceleration, the total load is linearly 
distributed along the keel and the center of pressure is located at a 
distance equal to one— third the wetted length forward of the step. 

Under actual three-dimensional conditions, however, as a result of 
the longitudinal components of flew introduced by the pressure gradient 
along the keel and the finite length, the theoretical two-dimensional 
distributions are somewhat modified as qualitatively shown in figure 2 
by the broken— line curves, so that 

F = [fO)] 2 [f 1 (A,s)z 2 z + 2f 2 (A,s)zz2] 

where the reduction in force from that calculated by the assumption of 
two-dimensional flow in the flow planes is determined by the geometry 
of the Immersed 'part of the hull and the station under consideration. 

As a result of this effect the total load on the seaplane was 
reduced by the application cf an aspect— ratio type of correction 0(A) 
which depends on the immersed shape of the hull. In the present analysis 
a similar type of reduction 0^(A) is applied to the moments calculated 
by assuming two-dimensional flow in Hie flow planes. The pitching 
moment about the step-keel point is therefore given by the expression 


Mb = [f(f>)] 2 j^(A) ^ J 


z s /tan t 


(z 2 ^ + 2zz 2 )(l — s) ds 


or 


^ - 


[f(p)J\(A) PI fak 


+ z 2 z b 3 


6 tan^T 


k 


( 12 ) 
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where tlie pitching moment is taken as positive In the direction of 
increasing trim. 

The pitching moment may he related to the motion of the seaplane 
measured relative to the water surface by introducing the following 
expressions from reference 2: 


where 


and 


2 s = 


7 


COS T 


Z = 


+ K 


COS T 


( 13 ) 

( 1 *) 


K = (xq — y Q tan t ) sin 


z = 


( 15 ) 


COS T 

Substituting equations (13), (l 4 ), and (15) into equation ( 12 ) gives 


M s = 



77 k 

6 sin^r oob^t 

_4 


+ Cy + K COS T ) 2 y3 


(16) 


Equation (l6) may be reduced to dimensionless form by introduction 
of the nondimens ional variables defined by equations (5), (6), (7), (8), 
and ( 9 ). Thus, the equation 



(IT) 


gives the relationship between the pitching-moment coefficient at any 
instant and the corresponding draft coefficient, vertical— velocity ratio, 
and load— factor coefficient, where the pitching-moment coefficient for 
the step-keel point is defined by 



Sin T cos T 


(18) 
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Equations (l^), (l 6 ), and (l 8 ) of reference 2 provide the following 
relationships which define the motion of the seaplane during the impact: 


C 


2 = 



1 + C 4 3 


(19) 



Combining equations ( 17 ), ( 19 ), and (20) results in the following equations 
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Applying equations (20 ) and (21 ) provides the relationship between 
the pitching-moment coefficient and the vertical-velocity ratio: 



The preceding equations apply at all instants during an impact. 


Center of Pressure 

The cente 2 >-of— pressure distance p from, the step is defined by 



(26) 
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From equations ( 5 ) and (l8), since 



W 7 

g COS T 



COS T 



and 

M ^ il (a) 

0 sin t cos t \g / 0(A) 


the relationship that exists at any instant between the dimensionless 
center— of -pres sure coefficient C C p and the moment and load-factor 
coefficients is given by 


where 



w 


( 27 ) 


'cp 


= p sin t 



0( A ) (r ]&(P)] 2 0(A)p7t| 
.01 (A.) 


^/3 


6 sin t cob2t 


( 28 ) 


The combination of . equation (22) with ( 27 ) gives the relationship 
between the center^-of— pressure coefficient and the draft coefficient: 


C C p 




(29) 
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Comb ining equations ( 2 ^-) and. (27) gives the foil awing r slat ionship 
between, the center- of— pressure coefficient and. the vertical-velocity 
ratio: 



Substituting the geometric relationship X = 


7 


sin t 


(see fig. 2) and 


the definition of the draft coefficient, equation (6), into equations ( 27 ) 
and (28) gives the following expression for the ratio of the center-of— 
pressure distance to the wetted length: 


( ^m s 

* = Wi 


(31) 


where 


p 0 (A) 


Combining equation (31) with ( 20 ), ( 21 ), and (25) provides the 
relationship between the ratio of the c ent er— of — pr a s sure distance to the 
wetted length and the vertical-velocity ratio which exists at all instants 
during the impact: 


C r - 


1 _l_ 1 + k 

^ 12 + K 

3^0 



( 3 &) 
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Combining equations (20 ) and (32) gives the relationship "between the 
ratio of the c enter— of -^pre e sure distance to the wetted length and the 
draft coefficient which applies for all values of ie: 

°r = 3 + ]| 0 a 3 ■ (33) 


Transfer of Moments 

The determination of the pitching moment and the center of pressure 
permits the ready calculation of the hydrodynamic moment about any point 
on the seaplane. For any point o located at a distance a forward of 
the step (see fig. l) 


= F t (p - a) 


Introducing equation ( 26 ) gives 


M~ = M a /1 - - 
°h ^ P/ 


By use of equations (l 8 ) and (28)^ equation (34) may be written 


or 


Cm 


°h 


Cm Q 

<r~ 

cp 


'cp 


a aln T L f a 0(A)p»t ^ 

t(A) [ 6 sin t cos 2 rj/ 


where 


P M °h /g\ 0 ( A ) 

^ = ;2 0l(A ) 


(34) 


(35) 


(36) 


sin t cos t 


(37) 
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Special Conditions 

Conditions at maximum moment .— Equation (25) provides the general 
relationship "between the pitching-moment coefficient and the vertical- 
velocity ratio which applies at all times during an impact. The condi- 
tions which exist at the instant when the maxi mum moment is reached 

d Cm 0 

are determined "by differentiating equation (25) and setting = 0. 

Thus, the equation 


1 + k 


_ 2 _ 

fo 


+ K 



+ 


2 (h + K )-^T^ K t + x6KS 


~ + 2k 
Yo 


0 


( 38 ) 


provides the relationship "between the vertical— velocity ratio at the 
instant of maximum moment and the approach parameter k . The substitution 

of values of obtained by solution of equation (38) into equations (2^), 

Yo 

(30), (20), and (2l) permits the determination of the maximum pitching- 
moment and center— of ^pressure coefficients as well as the draft and load- 
factor coefficients at the instant of maximum moment. The time coefficient 
corresponding to the occurrence of maximum moment may be calculated from 

the variation of C+ with determined in reference 2 . 

Yo 

Conditions at maximum acceleration .— The relationship between the 
vertical— velocity ratio at the instant of rnwylmn-m acceleration and the 
approach parameter tc is given by equation (27) of reference 2 : 


1 + JC 



- logo 




+ 6k 


(1 + k) 



( 39 ) 


The vertical-velocity ratio at the instant of maximum acceleration may- 
be determined from the solution of equation ( 39 ), and the maximum load- 
factor coefficient and the draft coefficient at maximum acceleration 
(reference 2 ) may be calculated by means of equations ( 2 l) and ( 20 ). 
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By use of these values of the vortical— velocity ratio, the pitchlng-sncsoent 
and. center— of— pressure coefficients at the instant of maximum, acceleration 
may he calculated from equations ( 25 ) and ( 30 ). 

y 

Conditions at maximum drafts At the instant of maximum draft — - 0 . 

y Q 

With this substitution equations (20) and (21 ) permit the determination of 

the maximum draft and of the load— factor coefficient corresponding to this 
instant (reference 2). Similarly the application of equations ( 25 ) and 
( 30 ) provides the following explicit relationships between the pitching- 
moment and center-of— pressure coefficients at the instant of maximum 
draft and the approach parameter k : 


Cme 



K 

1 + K 



(40) 


and 


C 


_ 1 

cp ~ 12 


f 

1 4- K 
* 


1 

1+ It 



(40a) 


The expression relating k and the ratio of the center-of— pressure 
distance to the wetted length at the instant of maximum draft is obtained 
from equation ( 32 ): 


C r = £ + 


1 1 + rt 

12 k 



(41) 


Limiting con d itions .— Since the approach parameter k may range 
between 0 and it is desirable to determine the limiting values 

between which the coefficients of motion at different stages of the impact 
may vary. The condition of k = 0 is obtained when the flight path at 
contact is normal to the keel of the seaplane. For this condition 
equations (19) and (20 ) become 



1 + C d 3 


( 42 ) 
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and 


i-i- 

c d 3 = r— a 


(43) 


Combining equations (42) and (43) gives 




m 


and 


G l = 


30j 

(l + °a 3 ) 3 


The time coefficient may "be directly obtained by integrating 
equation (43) I 


C t = C d (l + |c d 3) 


(45) 


(46) 


or 



y\i/3 



For the case -where K = 0, equation ( 25 ) becomes 


(47) 




m 


Combining equations (43) and (48) gives 

e a 3 (4 + °a 3 ) 


Cm s 


4(i + °a 3 )~ 


(49) 
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For this condition equation (30) gives 



The relationship between the center— of— pressure coefficient and the draft 
coefficient is given by equation ( 29 ) which is valid for a l l values of k. 

For k - 0 the variation of the ratio of the center— of -pressure 
distance to the wetted length is obtained from equation ( 32 ) 





(51) 


The preceding equations apply at all instants during an impact in 
which the resultant velocity is normal to the keel and permit the 
determination of the conditions at specific stages of the impact which 
may be of interest. For k = 0, equation (38) pezmiits the calculation 
of the vertical— velocity ratio corresponding to the occurrence of 
maximum, moment: 


7 


r o /13 - 2 


= 0.6228 


(52) 


The draft coefficient at this instant is determined from equation (43): 


c d = (1/I3 - 3) 1/3 = 0. 


8460 


(53) 


while the load— factor coefficient is obtained from equation (44): 


C 


l 


3(1/13 - 3) 2 / 3 

(tfj - 2) 3 


0.5187 


(5^) 
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Equation (47 ) permit s the calculation of the value of the time coefficient 
at -which the maximum moment occurs : 

C t = (jft3 - 3) 1//3 + |(i/l3 - 3)^ 3 = 0.9741 (55) 

The maximum moment coefficient is determined from equation (48): 



(1/I3 + l)(i/l3 - 3) 

4(|/5J - 2) 3 


0.1685 


(56) 


■while equation ( 50 ) permits the calculation of the center— of— pressure 
coefficient for the same instant: 




(]/l3 + l) = 0.32^7 


(57) 


The ratio of the center— of— pressure distance to the wetted length is 
obtained from equation ( 51 ) : 


C*. = 


i^3 h 1 

12 


0.3838 


( 58 ) 


In a similar manner 3 the solution of equation (39) permits the 
determination of the conditions at the instant of maximum acceleration. 
As shown by equations ( 38 ) to (4l) of reference 2 the state of motion at 
maximum acceleration for k = 0 is defined by the following constants: 


C 


d 


/o\l/ 3 

(f) - 0.6586 


(59) 


7_ 

fo 


7 

9 


( 60 ) 


iif%f - °- 6123 


( 61 ) 
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and 



0.7057 


( 62 ) 


Equation ( 60 ) penults the determination of the moment coefficient 
at maximum acceleration from equation (48). Eor k = 0 



0.1440 


(63) 


The center—of —pressure coeff icient^at this instant is obtained from 
equation ( 50 ): 



and the ratio of the center— of— pressure distance to the vetted length 
from equation ( 51 ): 


C r 



s: = °' 3571 


( 65 ) 


"When K = 0 equation (46) shows that the draft always increases 
with time. Thus, while the downward velocity of the float grows smaller 
as the impact progresses, a maximum, draft is never reached. As discussed 
in reference 2, this result is due to the fact' that, since the decelera- 
tion of the seaplane, for k - 0, is in the same direction as the 
resultant velocity (normal to the keel), the seaplane continues along its 
original path of motion throughout the impact. Thus, only the flow planes 
directly beneath the keel are affected by the Immersion and absorb all 
the momentum lost by the seaplane. Consequently, an infinite virtual 
mass (infinite draft) is required to satisfy the condition of zero vertical 
velocity of the seaplane. This result is due to the neglect of the 
buoyant forces which, because of the large drafts attained, are of 
greatest importance at the very high flight— path angles beyond the 
range of approach conditions applicable to conventional seaplanes. Eor 
values of rc other than zero the foregoing result does not apply since 
part of the momentum lost by the seaplane is contained in the downwash 
left behind the step, resulting in the attainment of a finite maximum 
draft (reference 2). 
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For small values of k (steep impacts) the differences "between the 
time coefficients corresponding to the occurrence of maximum acceleration, 
maximum moment, maximum draft, and the rehound from the water surface are 
large. As the flight path "becomes flatter this difference becomes smaller 
and approaches zero as the planing condition (k = co ) is reached. For 
the planing condition the dimensionless variables approach the following 
limiting values: 


lim 

Cl = 

K — >00 


lim 

y = 


Jo 

K — >00 


lim 

Cd - 

K — ^00 


lim 

Ct = 

K 


lim 

Croe = 

re — >00 


lim 

C cp = 


00 

1 

0 

0 


lim C_ = L 

r 3 

tc —^00 


THEORETICAL RESULTS AND DISCUSSION 


On the basis of the foregoing analysis the variations during an 
impact of the dimensionless coefficients representing the pitching 
moment and center of pressure depend solely on the magnitude of the 
approach parameter k. For a given value of re, therefore, the 
respective variations of the pitching-moment and center— of— pressure 
coefficients may each be represented by a single curve regardless of 
the seaplane properties, attitude,, or initial velocity. Consequently 
a single variation exists between the coefficients corresponding to any 

given stage of the impact |a given value of and the approach 

parameter k. The use of dimensionless variables, by thus talcing int o 
account such factors as dead— rise angle, weight, trim angle, and velocity 
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in accordance with, the laws governing the variation of the pitching 
moment and center of pressure with these quantities, permits reduction 
of the impact conditions to a common "basis defined "by the approach 
parameter. 

As shown "by the form of the coefficients, for a given value of k 
the pitching moment and center of pressure corresponding to a given stage 
of the impact are related to the primary variables constituting the 
float properties, attitude, and magnitude of the initial velocity "by the 
following proportionalities; 


M 


s 



sin t cos t \g/ 


£l(a) 

0(A) 


and is independent of the angle of -dead rise, and 


As a result of the analysis, figures 4 to 13 show the theoretical 
variations of the pitching-moment coefficient, centei^-of— pressure coeffi- 
cient, and ratio of the center*— of— pressure distance to the wetted length 
obtained during step— landing impacts. The solutions include values of tc 
which correspond to a wide range of approach conditions extending from 
impacts in which the resultant velocity is normal to the keel (tc = 0) to 
impacts along very shallow flight paths which begin to approach the 
planing condition. Figures l4 to 20 show the variations with the 
approach parameter of the conditions which exist at different stages of 
the impact i In order to compare the conditions at the instant of 
maximum pitching moment with those existing at other stages of the 
immersion, solutions for the state of motion and the time corresponding 
to the occurrence of maximum acceleration, maximum draft, and at the 
instant of zero draft during rebound, which were previously presented 
in reference 2, are shown in figures 17 to 20 by the broken— line curves. 

Although the dimensionless curves permit the complete determination 
of the pitching moments as well as the motion experienced by a seaplane 
during an impact, some interpretation of the results is desirable. For 
impacts at different values of k, since both the pitching-moment coeffi- 
cient and the time coefficient are based on the initial vertical velocity, 
the actual time histories of the moment will have the same relative 
shapes as the dimensionless curves if the vertical velocity is the same 
for each value of k . Thus the curves shown in figure 7 E^y he interpreted 


W — / ^ 

ffO)f0(A) p tan 2 T 0(A) 
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as corresponding to different seaplanes landing with, the same sinking 
speed "but with different resultant velocities and thus varying flight- 
path angles, in which case the maximum moment for the shallow approaches 
(high resultant velocities and large values of k) will he greater than 
that for the steeper flight— path angles and will he attained in a shorter 
time after contact. 

On the other hand, a somewhat different interpretation of the 
dimensionless curves may he given when a particular seaplane landing over 
a range of flight— path angles is considered. In this case the resultant 
velocity is more or less constant while the sinking speed, which depends 
largely on piloting technique, may he varied to provide a range of 
flight— path angles. As has been previously shown, the pitching moment 
at any proportional part of the impact cycle varies as the square of the 
initial vertical velocity while the corresponding time is inversely 
proportional, to the vertical velocity. If the resultant velocity and the 
trim angle are held constant, steeper flight paths are associated with 
the smaller values of k . Since the increase in vertical velocity more 
than offsets the reduction in pitching-moment coefficient and the increase 
in time coefficient with decreasing values of k, the maximum moment 
obtained with constant resultant velocity and trim angle will he greater 
for the steeper approaches than at the low flight— path angles and will 
he attained in a shorter time after contact. The increase in pitching 
moment with flight-path angle is illustrated by figure 21 which shows the 
variation with flight— path angle and trim of a dimensionless pitching- 
moment coefficient based on the resultant velocity at contact. 

With regard to the sequence of events during an impact, figures l4 
to 20 show that the maximum moment occurs slightly after the maximum 
acceleration is reached and precedes the attainment of the maximum draft. 
The state of motion at the Instant of ma xi mum moment Is only slightly 
different from the conditions which exist at the instant of maximum 
acceleration. As might he reasonably expected, for the same sinking 
speed, a greater time is required to reach a given stage of the impact 
at the high flight— path angles than is required for the flatter approach 
conditions. In a similar maimer, the differences in the states of motion 
and the times corresponding to the various stages of the impact are large 
for small values of k (high flight— path angles) and become very ema 1 1 
as the planing condition is approached. 

Although the pitching moment at any stage of an impact varies with 
the square of the initial vertical velocity, the location of the center 
of pressure is Independent of the initial velocity and Is determined 
primarily by the wetted length. In fact, as is shown by equation (33)j 
the centeD>-of— pressure distance forward of the step is only slightly 
greater than one— third the wetted length. The fact that this distance Is 
somewhat greater than one— third the wetted length arises from the 
quadratic nature of the longitudinal distribution of the negative 
increment in hydrodynamic load accompanying the deceleration of the 
Yirtual mass. (See fig, 2.) When this distribution Is added vectorially 
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to the positive linear load variation caused by the expansion of the 
virtual mass, the resulting distribution of the total load is not quite 
linear and the center of pressure is shifted slightly forward of the 
center— of— pressure location for a linear distribution. The extent of 
this forward shift of the center of pressure depends, of course, on 
the relative magnitudes of- the quadratic and linear distributions and 
increases with penetration (see equation (33))* Consequently, th.6 most- 
forward location of the center of pressure, both in an absolute sense 
and in terms of the wetted length, for any given value of k will be 
attained at the instant -of maximum draft. For impacts at low values 
of k this forward shift of the center of pressure in terms of the 
wetted length may be quite large as a result of the large drafts and 
the relative magnitudes of the quadratic and linear distributions 
obtained at the high flight-path angles. For the planing condition, 
since there is no acceleration, the total load is linearly distributed 
along the keel and the center of pressure is located at a distance equal 
to exactly one— third the wetted keel length forward of the step. As 
shown by figures 13 and 1 6 , even for the practical range of seaplane 
impact conditions (values of k > 0.2) the shift 1 of the center of pressure 
forward of the one— third point is relatively small and for practical 
purposes may be neglected in most caseB. 


COMPARISON WITH EXPERIMENTAL DATA 


The applicability of the theoretical results is illustrated by 
comparisons with experimental data obtained in the Langley impact basin, 
under controlled conditions in smooth water, with two float forebodies 
having angles of dead rise of 30° and 40° at the step. The total 
weights in the tests ranged between opproximately 1200 and 1350 pounds. 

The detailed test conditions are presented in table I. During the tests 
wing lift was simulated by the action of a pneumatic cylinder and cam 
system which was designed to apply a constant upward force to the float 
equal to the total weight. 

The experimental results include measurements of the horizontal and 
vertical components of velocity as well as time histories of the pitching 
moment. The moments were measured by means of the strainr-gage dynamometer 
truss schematically illustrated in figure 1 and are referred to the front 
float attachment point as shown in the figure. The original oscillograph 
records from a typical test-are reproduced, greatly reduced, in figure 22. 

Experimental pitching-moment- time histories, as derived from the 
oscillograph records, are shown in figures 23 to 2 6 . These data are 
compared with theoretical time histories of the total pitching moment 
about the front attachment point as well as with theoretical time 
histories of the hydrodynamic pitching moment about this point. 

The time histories of the hydrodynamic pitching moment were 
calculated by application of equations (35) and (37) in conjunction 
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with tlie analytical results plotted in figures 7 and 12. The dimension a 
in equation (35)* 'which represents the distance from, the step to the 
front attachment point, was measured as 2.89 feet for the two floats 
tested. In converting from, dimensionless variables to the dimensional 
quantities moment and time, it is necessary to introduce values for the 
functions f(p), 0 (a), and 0 i(a). Equations (10) and (ll) present 
expressions for- f(p) and 0 (a) which have been shown to be valid 
for a wide range of dead-rise and trim angles. Since 0 ^(a) has not 

been evaluated at the present time, it was assumed that the aspect— ratio 
correction to the total load 0 (a) applies uniformly to all flow 
planes; that is,- 0 j_(a) = 0(a). From the nature of this approximation 

it would appear that this assumption does not involve serious errors. 

In order to properly compare the theoretical and experimental 
results, the inertia and static moments, introduced by the fact that 
the center of gravity of the model did not coincide with the center of 
moments, must be added to the hydrodynamic moment about the front 
attachment point. The increment in pitching moment arising from this 
source is given by the equation 


Mo 


= W.p £ 

L 6 


COS T 


— (b cos t + c sin t) 


( 66 ) 


Thus the total pitching moment referred to the point o is expressed by 


Mot = 


g COS T 



1/3 


C^W 01 (A) e g |?(Pl 0(A)pir [\ , \ 

sin t 0(a) ^ l W i o_ I / V f/ 


6 sin t cos^r 


— VT ; p(b cos t + c sin t) ( 67 ) 

where the dimensions b and c define the location of the center of 
gravity relative to the center of moments, as shown in figure 1. Values 
of b and’ c for each of the configurations tested are given in 
table I. 

Figures 23 (a), (b), and (c) show the results of three tests of the 
float with 30° angle of dead rise at 12° trim. All three tests were 
made at approximately the same flight— path angle of about 0 and, 
therefore, closely correspond to a single value of the approach parameter k 
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equal to about 2.2. Two of the 'teste were made with almost identical 
resultant velocities. The third test was made at a velocity approxi- 
mately 60 percent greater than that for the first two runs. The 
theoretical results were calculated- for values of the initial vertical 
velocity measured in the teste and for a value of k = 2 , which 
approximately corresponds to the magnitude of the approach parameter 
associated with the tests. The theoretical total pitching moments 
appear to he in substantial agreement with the results measured in th.9 
tests. 

Figure 2k is a composite plot of the quantity 


+ Wf(b cos t + c sin r) 

_ 

g 


Cm g 0j(A) 

Bln T 0(A) 




[f(p)] a 0(A)p>t 


a/3 


sin t 


cos^t | 



(67a) 


against the variable ty Q . For the three tests previously discussed, 
since all parameters except the resultant velocity were held constant, 
equation ( 67 a} and the form of the time coefficient show that the results 
of all three tests should be reduced to the single theoretical variation 
given by the solid— line curve. As is evident from th .9 figure, this 
result is closely attained. The extent of the deviations which do exist 
may be taken as a direct indication of the consistency of the experi- 
mental data. The oscillatory nature of the results of some of the tests 
(see fig. 22 ) is attributed to structural vibrations induced in the 
equipment by the magnitude of the catapulting accelerations required 
to produce the relatively high horizontal velocities attained in the tests . 

Figures 25 and 2 6 show comparisons of theoretical and experimental 
pitching-moment time histories for typical impacts of a float with a 
40° angle of dead rise at trim angles of 9° and 6 ° and values of the 
approach parameter ic approximately equal to 1. Comparisons between 
theoretical and experimental maxLmum total pitching moments for this 
float are presented in figures 27 (a), (b), and (c) for trim angles 
of 12°, 9°, and 6 °. In order to reduce the increment in pitching 
moment due to the displacement of the center of gravity from the center 
of moments, another , set of tests was made, at 12 ° trim, with the float 
weighted to move the center of gravity rearward to a position vertically 
in line with the front attachment fitting. The results of these tests 
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are shown In figure 28. As is evident from, the figures, the theoretical 
results appear to "be in good agreement with the experimental data. 

A more exact evaluation of the end— flow correction to the pitching 
moment 0p(A), which in this comparison has "been taken equal to the 
end— flow correction to the total load 0(A), should result in even tetter 
agreement between the calculated and experimental results. In the absence 
of suitable pressure— distribution data, the favorable results of this 
comparison ma y be considered an indirect indication of the validity of 
the longitudinal distribution of the hydrodynamic load (running load) 
specified by the theory. 


SUMMARY OF EESUIZTS 
Theoretical 


The preceding analysis of the pitching moments experienced by 
Y— bottcm seaplanes during step— landing impacts has shown that: 

1. The pitching moment about the step and the cent er-of— pressure 
distance forward of the step may be represented in generalized form by 
means of the following dimensionless variables: pitching-moment 

coefficient 


„ M a /g\ ?(A) 

^ ‘ j 0 Z W fcCA) 

and center— of— pressure coefficient 


sin t cos t 


-cp 


= p sin t 


g J [f(p)l 2 0(A)pJ V 


6 sin 


It is also shown that the variation of these dimensionless quantities 
during an impact is governed solely by the magnitude of the approach 
parameter 


K = 


sin r 
sin 7 ,- 


cos 


(t 


?o> 


which depends only on the trim and the flight path at the instant of 
initial contact with the water surface. For a given value of k, 
therefore, the respective variations of the pitching-moment and center— 
of— pressure coefficients may each be represented by a single curve 
regardless of what the seaplane properties, attitude, or initial condi- 
tions may be. Consequently a single variation exists between the 
coefficients existing at any stage of the impact and the approach 
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parameter k. The expression of the pitching moment and center of 
pressure in terms of dimensionless variables, by thus taking into account 
such factors as angle of dead rise, weight, trim angle, and initial 
velocity, permits reduction of all impact conditions to a common basis 
defined by the approach parameter and leads to a simplified presentation 
of both theoretical and experimental results. 

2 . Item the form of the pitching-moment coefficient it can be seen 
that the pitching moment about the step at any given stage of the impact 
is independent of the angle of dead rise. 

3 . The maximum pitching moment about the step is attained slightly 
after the maximum acceleration (load) is reached and occurs prior to the 
attainment of the maximum draft. The state of motion of the seaplane 
and the time corresponding to the instant of maximum pitching moment are 
only slightly different from the conditions which exist at the instant 

of maximum acceleration. The differences in the coefficients representing 
the pitching moments, the states of motion, and the times corresponding 

to given stages of the impact (given values of such as the instants 

\ To/ 

of maximum acceleration, maximum pitching moment, maximum draft, or zero 
draft 'during rebound, are greatest at the low values of the approach 
parameter k (steep flight paths for conventional trim angles) and 
decrease as the .planing condition is approached. 

4 . For given seaplane properties, attitude, and constant vertical 
velocity, the pitching moments about the step are greater for the higher 
values of the approach parameter k (corresponding to low flight— path 
angles and high resultant velocities, for conventional trim angles ) than 
those at the lower values of k (steep flight paths and low resultant 
velocities, for conventional trim angles) and are attained in a' shorter 
time after contact. For a given resultant velocity, on the other hand, 
the greater pitching moments are attained at the lower values of K . 

5. In a step Impact the center of pressure is located at a distance 
only slightly greater than one-third the 'wetted keel length forward of 
the step. ' The relationship between the cent ei^-of-pres sure distance and 
the wetted length is independent of the magnitude of the initial velocity 
and is determined pr im a r ily by the wetted length. For tb .9 planing condi- 
tion the center of pressure is located at a distance exactly equal to 
one— third the wetted keel length forward of the step and is shifted 
slightly forward with decreasing values of the approach parameter k . 

The most forward center— of ^pressure locations, both in an absolute sense 
and in terms of the wetted length, are obtained in conjunction with the 
large draft coefficients reached at the low values of k (steep flight 
paths for conventional trim angles). For the practical range of impact 
conditions (k ^ 0 . 2 ), however, only small errors will be introduced if 
the resultant force is assumed to act at a distance equal to one— third 
the wetted keel length forward of the step. 
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Experimental 

Comparison of the theoretical results -with test data obtained in. 
the Langley Impact basin with, floats haying angles of dead rise of 30° 
and 40° at the step indicates that, within the limits of experimental 
error, the calculated pitching moments are in good agreement with th© 
experimental resulte * 

CONCLUDING- REMAKES 


A theoretical investigation has been made of the hydrodynamic 
pitching moments experienced "by Y— bottom seaplanes during step— landing 
impacts. Hie analysis shows that the pitching mom e nt and the center— 
of— pressure location, may he represented in generalized form by means of 
dimensionless variables, respectively designated the pitching-moment 
coefficient and th© center— of —pres sure coefficient, which- take into 
account such factors as th© seaplane wotght, angle of dead ripe, trim 
angle, and initial velocity. The variation of these coefficients 
during an impact is governed by the magnitude of the approach parameter 
determined by the trim and the initial flight— path angle. 

Equations are presented from which the variation of the pitching- 
moment and center^-of— pressure coefficients hays b<?ep calculated for a 
wide range of conditions extending frcm impacts along shallow flight- 
path angles approaching planing to very steep impacts in which the 
resultant yelogfty is normal to the keel. Solutions §r© also presented 
for the conditions which exist at th© instants of mxiTmim acceleration, 
Trwv-f rmi-rn pitching moment, maxlrnirm draft, and at the instant of zero draft 
during rebound, 

Comparisons of theoretic al pitching-momant time histories ash 
values of the maximum pitching moment with experimental data obtained 
in the Langley impact bPSip With floats of 3Q° and ^fO 6 angle of dead 
rise indicate that the calculat§h results jape in good agreement with the 
measured values. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Eield, Ya. , January 30^ 19^8 
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Figure l- 


Float- forebody model attached to boom of impact- basin 
carriage by means of dynamometer truss. 
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Figure 4* — Theoretical variation of pitching rnoment with load factor „ 
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f^iqure S . — Theoretical narration of pitching moment with vertical velocity. 
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Figure 6 . — Theoretical rariation af pitching moment with draft. 
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Figure 7. — Theoretical ra nation of pitching moment with time. 
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Figure 6. — Theoretical ranation of pitching moment mth center of pressure. 
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Figure 3. — Theoretical rarigtion of center of pressure with load factor. 
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Figure 1 1 . — 'Theoretical rar/artion of center of pressure with draft. 
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Figure 42 . — Theoretical ra nation of center of pressure with time . 




Pa-fio of cen+er-of-pnessor& 
to wei fed length , Cp ■= 


* 

,3 


.2 


./ 


O J .2 .3 .& 



Draft coefficient y Crf ~y 

Figure (3. — — Theoretical ra nation of ratio of 

to wetted length with draft. 
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Figure 14 . — Theoretical variation of pitching moment at 


various stages of fhe impact mfh approach parameter. 
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Figure /J. — Theoretical variation of center of pressure 
at various stages of the impact with approach 
parameter . 
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Figure /6.— Theoretical variation of ratio of cepter-of -pressure distance 
to wetted length at various stages of the impact with approach 
parameter . <o 
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Figure fl. Theoretical ra nation of load factor at 

various stages of the impact with approach parameter. 
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Figure 16 . — Theoretical ra nation of vertical velocity at various stages of the impact 

with approach para mete r+ 
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of the impact with approach parameter. 
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Figure 20 . — Theoretical /a nation of time at various 
stages of the impact with approach parameter 











Pitching moment f I Id - ft 


56 


NACA TN No. I 63 O 



t 

I 

■Q 




O' 

c 


(a) Experimental data t fy Q = d.&d feet per second $ 
Vfrj- S4.CS feet per second f '1 0 - S./S f C-2..20- 



77 me after contact t sec 



(id) Experimental data: tr a - C.2.6 feet per second $ 
Yh 0 ~ E7.S3 feet per second / 7 0 = .£,3/ /C= £.!$. 


Ei pare S3 . Comparison between theoretical and experimental time histones 

of pitching moment for a P- bottom seaplane with an anqte of dead 
rise of 30°. W- f£3! pounds ; P- f£°. 
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Figure 23 . Concluded. 



Figure 24 Comparison of theoretical and experimental variations of total 

pitching moment for a V~bottom seaplane with an angle of dead rise of 
30 0 and different initial conditions . W a 12.31 pounds j T*- 12°. 
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Figure 25 . -Comparison between theoretical and experimental time histones 
of pitching moment for a V-bottom seaplane with an angle of dead rise 
of 40°. W m /2/3 pounds ; r • 6' s V\i a * 8.67 feet per second 1 Vj )a * 85.47 feet 
per second ; 7 a » 5. 79 " j C a /.0 . 
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Figure 26.- Comparison between theoretical and experimental time histones 
of pitching moment ‘ for a V-bottom seaplane ' with an angle of dead rise 
of 40°. W ■ 1213 pounds j 7*9° ■ j Vy 0 = 6.75 feet per second ; Vp 0 * 63. 69 feet 
per second ) 4 0 = 1. 6? ° j c* I.IO . 
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Figure 2.7.- Comparison between theoretical and experimenta / 
variation of maximum total pitching moment w/fh 
approach parameter. / 3 * 40 °; 1213 pounds > Wf =3Sopounds t 

a= 2.89 feet; 6 = 1.05 feet ; c= 1.8/ feet. 
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Figure za. - Continued. 
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Figure 27 .- Cone I tided . 
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Figure 28 . - Comparison between theoretical^ and experi- 
mental variation of maximum total pitching moment 
with approach parameter . /3=*o° s 12 0 

W~/343 pounds i Wf ~ 590 pounds a= 2.89 feet ; 
b°-.28feet • c* 1.4 6 feet. 



